Thermal Expansion Coefficients of Bi2Se3 and Sb2Te3 Crystals from 10 K to 270 K 
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Lattice constant of Bi2Se3 and Sb2Te3 crystals is determined by X-ray powder diffrac- 
tion measurement in a wide temperature range. Linear thermal expansion coefficients 
(a) of the crystals are extracted, and considerable anisotropy between a\\ and a± is 
observed. The low temperature values of a can be fit well by the Debye model, 
while an anomalous behavior at above 150 K is evidenced and explained. Griineisen 
parameters of the materials are also estimated at room temperature. 
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Recently, much attention has been given to an intriguing class of materials, the so-called 
topological insulators (TIs). This type of material exhibits a band gap in the bulk, but 
gapless states on the edge or surface, which are protected by topological order and cannot 
be analogized to conventional semiconductors or insulators^'^. Bi2Se3, Bi2Te3 and Sb2Te3 
are among the most interested compounds of three-dimensional TIs, owing to their robust 
and simple surface states^. Although these compounds were under extensive studies in 1950s 
and 1960s as excellent thermoelectric materials, some basic physical properties still remain 
unexplored. In this letter, we present the measurements of the temperature dependent linear 
thermal expansion coefficients of Bi2Se3 and Sb2Te3 crystals using X-ray powder diffraction 
(XRD). Thermal expansion is the tendency of materials to change in size and shape as they 
heat and cool. It is essential to device design and engineering, as the induced strain could 
cause the deformation of the device and affect its phonon dynamics. Indeed, our recent 
Raman spectroscopy study of TIs has uncovered significant contributions in the tempera- 
ture dependent phonon frequency shifts from the thermal expansion of the materials^. In 
addition, the knowledge of thermal expansion coefficients is necessary for the directional 
growth of TI crystals and the understanding of the high thermoelectric efficiency^. 

Large grain polycrystalline Bi2Se3 materials (single crystal grain size >1 mm) were syn- 
thesized at Sandia National Laboratories. First, Bi2Se3 pieces (99.999%, from VWR inter- 
national, LLC.) were placed in an evacuated (<10~^ Torr) quartz ampoule and melted at 
800 °C for 16 hours. The melt was then cooled at 10 °C/h to 550 °C, held for 3 days at 
this temperature, and finally allowed to furnace cool to room temperature. Single crystals 
of Sb2Te3 were grown by Bridgman method at Purdue University. Stoichiometric amount 
of high purity antimony and tellurium (99.999%) was deoxidized and purified by multiple 
vacuum distillations under dynamic vacuum of <10^^ Torr, and then heated up to 900 °C. 
This was followed by a slow cool down under a controlled pressure to minimize tellurium 
defects. Afterwards, the crystals were grown at a speed of 0.5-1.5 mm/h with a linear tem- 
perature gradient set to 5 °C/cm. Bi2Se3 and Sb2Te3 crystals have similar rhombohedral 
structure with five atoms in the trigonal primitive cell. A straightforward way to visualize 
the structure is to use a hexagonal lattice with the unit cell being a quintuple layer, as shown 
in Fig. 1(a). Like graphite, adjacent Se-Se (Te-Te) layers are hold together by weak van der 
Waals force. 

The XRD patterns were recorded using a Ruber G670 imaging-plate Guinier camera 
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FIG. 1. (color online) (a) Hexagonal unit cell of Bi2Se3. The distance between two nearest Se 
atoms of the two adjacent Se layers is larger than the sum of their covalent radii. As a result, the 
Se-Bi-Se-Bi-Se quintuple layers are hold together by weak van der Waals force. Bi2Te3 and Sb2Te3 
have a similar structure by substituting Se atoms with Te, Bi atoms with Sb. Right panel: XRD 
spectra of (b) Bi2Se3 and (c) Sb2Te3 at 10 K. Maximum peak intensity is normalized to 1. 



equipped with a Ge monochromator and Cu K^i radiation (1.54059 A). Data were collected 
in steps of 0.005° in a wide temperature range from 10 K to 270 K. The lattice parameters 
were calculated via WINPREP program with residual factor 0.02 < < 0.03. Typi- 
cal XRD spectra of Bi2Se3 and Sb2Te3 crystals (at 10 K) are shown in Fig. 1(b) and (c), 
respectively. From the spectra, one can determine the lattice constant of Bi2Se3 (extrap- 
olated to K): Ohcx = 4.1263 A and Chcx = 28.481 A, which translate to aj-ho = 9.7880 A 
and arho = 24.337° for the rhombohedral cell. For Sb2Te3, the corresponding values are 
Obex = 4.2423 A, Chcx = 30.191 A, and a^ho = 10.357 A, a,ho = 23.635°. 

Figure 2 plots the measured lattice parameters of Bi2Se3 and Sb2Te3 as a function of 
temperature, from which the linear thermal expansion coefficients (a) of the material can 
be derived 

1 

(1) 



aL 



1 dL 



where L is the lattice constant (i.e., ahex or Chex), and T is the temperature. Figure 3 shows 
the obtained linear thermal expansion coefficients of Bi2Se3 and Sb2Te3 parallel (an) and 
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FIG. 2. Lattice constant of (a) Bi2Se3 and (b) Sb2Te3 as a function of temperature. 

perpendicular (aj,) to the trigonal c-axis. Considerable anisotropy between a\\ and a±_ is 
observed, reflective of the material bonding anisotropy. The ratio of the elastic constants 
C11/C13 characterizes the anisotropy of chemical bonds; Cn/Ci^ = 2.5 for Bi2Te3^ and 
similar values are expected for Bi2Se3 and Sb2Te3 owing to the structural similarity and 
comparable Debye temperature^. To understand the temperature dependence of a(T), it is 
instructive to introduce the Griineisen relation* 



Ci,7 = 



(2) 



where Cy is the specific heat, 7 is the Griineisen parameter, k is the bulk modulus, and 
P = 2a|j + aj_ is the volumetric thermal expansion coefficient. Usually, 7 and n are only 
weakly dependent on the temperature^, indicating that the thermal expansion coefficient has 
the same temperature dependence as the specific heat. This relationship (a(T) vs. c^(T)) 
remains valid for anisotropic materials^°, therefore we can apply Debye model to fit our data 



a = ao{—) / — 

Qd Jo (e^ - 1)^ 



(3) 



where ao is a temperature independent fitting parameter, and O/j is the Debye temperature. 
Good agreement is achieved in the low temperature regime, as shown by the solid lines in 
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FIG. 3. (color online) Linear thermal expansion coefficients of (a) Bi2Se3 and (b) Sb2Te3 as a 
function of temperature. Solid lines represent best fits using Debye model. The data well agree 
with the Debye law at low temperatures, while considerable deviation is evidenced at above 
150 K. 

Fig. 3, consistent with that observed in Bi2Te3^^. The obtained Debye temperature (from 
fitting both a\\ and a±) is Qo = 160 K for Bi2Se3 and Qd = 200 K for Sb2Te3, close to the 
hterature values (182 K and 160 K, respectively)^. At above 150 K, the experimental results 
deviate from the Debye model. Similar anomalous behavior is also evident in Bi2Te3^^, but 
the origin is not conclusive^^ . There are two possible explanations: (1) This is due to higher- 
order anharmonic effects. Generally, the anharmonicity of the bonding forces in a crystal 
is characterized by the Griineisen parameter. Higher-order anharmonic effects may lead 
to a nontrivial temperature dependence 7(T), giving rise to anomalous <y\\{T) and a±{T) 
through the Griineisen relation. This scenario is used to explain the thermal expansion of 
tellurium^"^'^^ and Bi2Te3^. (2) The deviation could be due to the breaking of the van der 
Waals bond between two Se-Se (Te-Te) layers, as proposed in Ref.[5] and [12] for Bi2Te3. The 
van der Waals binding energy is typically in the range of 0.4-4 kJ/mole, while the energy of 
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TABLE I. Thermal expansion coefficients and Griineisen parameters of Bi2Se3 and Sb2Te3 (at 270 



K), and Bi2Te3 (at room 


temperature) . 
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thermal motion at room temperature is about 2.5 kJ/mole, suggesting that the degradation 
of the van der Waals bond may cause the anomaly in a{T) at levitated temperatures. These 
two reasons may also explain the differences in a(T) between Bi2Se3 and Sb2Te3, but further 
experimental and theoretical work is needed to fully understand it. 

Finally, we summarize the room temperature thermal expansion parameters of Bi2Se3, 
Sb2Te3 and Bi2Te3 in Table I, as a reference for future study of TIs and device engineering. 
The thermal expansion coefficients for Bi2Se3 and Sb2Te3 are taken from Fig.3, which are 
in good agreement with the values documented in the materials handbook^. Griineisen 
parameters are calculated from Eq. [2] . Specific heat at constant pressure (cp) of the materials 
is better documented in the literature, and thus used in the calculation^^. For Bi2Te3 the 
values are extracted from various previous works. 

In conclusion, we have measured the thermal expansion parameters of Bi2Se3 and Sb2Te3 
crystals in a wide temperature range from 10 K to 270 K. The extracted linear thermal 
expansion coefficients are found in consistent with the Debye model at low temperatures, 
but deviate at above 150 K. Our result is crucial for interpreting the temperature dependent 
Raman shift in Bi2Se3 and Sb2Te3^, and may provide insight for understanding the thermal 
properties of TIs. 
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